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We describe a recently developed generalisation of the Poincare sphere method, to represent 
pure states of a three-level quantum system in a convenient geometrical manner. The construction 
depends on the properties of the group SU(3) and its generators in the defining representation, 
and uses geometrical objects and operations in an eight dimensional real Euclidean space. This 
construction is then used to develop a generalisation of the well known Pancharatnam geometric 
phase formula, for evolution of a three-level system along a geodesic triangle in state space. 
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Sometime after the discovery by Berry in 1984 of the quantum mechanical geometric phase in the framework of 
cyclic adiabatic evolution [Q, Ramaseshan and Nityananda pointed out in an important paper Hthat as early as 
1956 Pancharatnam had put forward closely related ideas in the context of polarisation optics |J. Subsequently 
Berry himself analysed the significance of Pancharatnam's work in the light of later developments ||. The more 
recent quantum kinematic approach to the geometric phase brings out in fully transparent fashion the way in which 
a phase essentially identified by Pancharatnam is one of the two basic ingredients involved in the very definition of 
the quantum geometric phase, the other being the so called dynamical phase ||. 



a: 



Pancharatnam's work made essential use of the Poincare sphere representation for the manifold of pure polarisation 
states of a plane electromagnetic wave |1. As is well known, diametrically opposite points on the Poincare sphere 



correspond to mutually orthogonal polarisation states incapable of interfering with one another. For two states of 
polarisation not mutually orthogonal in this sense Pancharatnam introduced a physically motivated convention or rule 
which would tell us when these two states are " in phase " , i.e., capable of interfering constructively to the maximum 
possible extent. More precisely, this relationship is defined at the level of field amplitudes mapping onto given points 

on the Poincare sphere. He then went on to show that this relation of being " in phase " is not transitive. That is, if 

we take three polarisation states A,B,C on the Poincare sphere, and arrange that the fields mapping onto A and B 
are "in phase" , and similarly those mapping onto B and C are "in phase" , then in general the fields mapping onto A 
and C are not "in phase" . He also calculated the extent to which these last two fields are "out of phase" and showed 

that this "phase difference" equals one half the solid angle on the Poincare sphere subtended by the spherical triangle 
ABC obtained by joining the vertices A, B and C by great circle arcs (geodesic arcs) on the sphere. 

This fundamental and early result of Pancharatnam has found a natural interpretation in the context of two-level 
quantum systems, for which the space of pure state density matrices is again the sphere S 2 . In the modern terminology 
for geometric phases, one is (most often) interested in cyclic evolution in the state space, and the calculation of the 
associated geometric phase. Evolution along a great circle arc on S 2 is particularly simple in that it can be generated 
by a constant (i.e., time- independent) Hamiltonian in such a way that the dynamical phase vanishes. A geodesic 
triangle on S 2 is then the simplest and most elementary yet nontrivial cyclic evolution one can imagine for a two-level 
system; it can be produced by a piecewise constant Hamiltonian, and the dynamical phase can be arranged to vanish 
identically. Then the meaning of Pancharatnam's result is that the resulting geometric phase is one half the solid 
angle on S 2 subtended by the triangle M . 
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We may note in passing that while Pancharatnam's original result refers principally to the vertices A, B, C of a 
spherical triangle on S 2 , and calculates the degree of nontransitiveness of the relation of two field amplitudes being 
"in phase" , in casting it into the modern geometric phase language we are equally concerned with the great circle 
arcs connecting these vertices, since we deal with continuous cyclic unitary or Hamiltonian evolution of pure quantum 
states. It is for such evolutions that geometric phases arc customarily calculated. Two more remarks are not out 
of place at this point. One is that in fact geometric phases can be perfectly well defined for noncyclic (and even 
nonunitary) evolutions, though in this paper we shall not be concerned with them ||. The other is that for piecewise 
geodesic and overall cyclic evolutions in a general quantum system, geometric phases get related to certain invariants 
introduced by Bargmann long ago J9|, and these invariants are in conception very close to Pancharatnam's original 
motivations. 

We shall refer to the "half the solid angle" result as the Pancharatnam formula for geometric phases for two-level 
systems undergoing piecewise geodesic cyclic evolution along a spherical triangle on S 2 . The main purposes of this 

paper are to: (i) present a recently developed generalisation of the Poincare sphere representation for pure states 
of three-level quantum systems JTCfl ; and (ii) to then obtain a generalisation of the Pancharatnam formula for such 
systems. 

The material of this paper is arranged as follows: Section II outlines the generalisation of the Poincare sphere for 

three-level systems. The Poincare sphere S 2 gets replaced by a certain four-dimensional simply connected region 
O contained wholly within the unit sphere S 7 in eight dimensional real Euclidean space. The transitive action of 
SU(3) on O, via eight dimensional orthogonal rotations, some intrinsic properties of 0, and a local coordinate system 
for O, are described so as to assist in forming a picture of this object. Section III recalls very briefly the main 
features of the quantum kinematic approach to the geometric phase of a general quantum system. Both the roles of 
ray space geodesies and the Bargmann invariants are highlighted. The case of two-level systems, and the statement 
of the Pancharatnam formula, are then given. It is pointed out that it is a fortunate circumstance that ray space 
geodesies and geodesies on S 2 happen to coincide in the case of two-level systems. Finally the general formula for 
the geometric phase for any cyclic evolution of the three-level system is given. Section IV discusses the properties of 
ray space geodesies for three-level systems, and their representation as curves in O. We find that the latter, while 
they are plane curves, are not geodesies in the geometrical sense on S 7 . We also find that it is possible to construct 
constant Hamiltonians which would give rise to evolution along any such geodesies. Section V puts together the 
ingredients of the previous Sections to develop the generalised Pancharatnam formula for three-level systems. This 
involves describing the most general geodesic triangle for a three-level system state space, and then computing its 
geometric phase. Whereas a geodesic triangle on S 2 involves three independent intrinsic parameters, for a geodesic 
triangle in O it turns out that four independent intrinsic parameters are needed. The geometric phase then depends 
on all four of these parameters, and this is borne out by the explicit formula for the phase. Section VI contains some 
concluding remarks. 



II. GENERALISATION OF THE POINCARE SPHERE REPRESENTATION FOR THREE-LEVEL 

SYSTEMS 

We recall very briefly the salient features of the Poincare sphere representation for two level systems, using through- 
out the notations and terminology of quantum mechanics JTTJ . We deal with a two-dimensional complex Hilbert space 
, unit vectors in which are denoted by ip, if>' ■ ■ ■. The density matrix corresponding to a pure state ip is given by 

the projection p = ipip^. Its expansion in terms of the Pauli matrices <7j leads to the Poincare sphere construction: 

p = H^ = -(1 + n.tr), 

p f = p 2 = p > 0, Tip = 1 & 

n* = n, n.n = 1 <S> n e S 2 . (2.1) 

Thus each pure state in the quantum mechanical sense, or normalised ray, corresponds in a one-to-one manner to a 
point on the 2-dimensional unit sphere S 2 embedded in Euclidean three dimensional space 1Z 3 . Since 

P = = ~ (1 + n.tr) , p' = ip'ip'^ = - (1 + n'.cr) 

Tr(p'p) = |(^»| 2 -i(l + n'.n) (2.2) 
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we see that diametrically opposite points on 5 2 correspond to mutually orthogonal rays or Hilbert space vectors. Here 
(if)',ip) is the inner product in TL^ 2 \ Finally if a vector ip G HS 2 ^ is subjected to a transformation u G SU(2), the 
representative point n€S 2 undergoes an orthogonal rotation belonging to 50(3): 

i>' =utp, ue SU{2) => 
n'j = R jk (u)nk, 

Rjk(u) = ^((TjWfeut), R(u) G SO{3). (2.3) 

As is well known, all elements 1Z G 50(3) are realised in this way and we have the coset space identifications 
S 2 = SU(2)/U(l) = 50(3)/50(2). 

Now we present the natural generalisation of this construction to three-level systems. We deal with a 3-dimensional 
complex Hilbert space TiS 3 ' , elements of which will be again denoted by ip,^', • ■ •■ The roles of SU(2) and the Pauli 
matrices Oj are now played by the group SU{3) via its defining representation, and the eight hcrmitian generators A r 
in this representation : 

SU(3) = {A = 3 x 3complex matrix \A^A = 1, det A = 1} ; (2.4a) 



Ai 
















, A 3 = 




-1 


• 


A 4 = ( 




: 











A 7 = -i , A 8 = —= 1 . (2.4b) 




The matrices A r obey characteristic commutation and anticommutation relations: 

4 

[A r , A s ] = 2if rst \t, {A r , X s } — g^rs + 2d rst \t] 

V3 1 

/l23 = 1, /458 = /678 = "77", /l47 = ^246 = /257 = ^345 = /5I6 = /637 = ^ > 
dns = ^228 = ^338 = — <^888 = _ Txi ^448 = <^558 = <^668 = <^778 



1*440 — "000 — — "-((O 2V3' 

dl46 = rfl57 = — ^247 = ^256 = <^344 = <^355 = — <*366 = ~ <^377 = (2-5) 

Here we have given the independent nonvanishing components of the completely antisymmetric f rs t and the completely 
symmetric d rst ; the former are the SU(3) structure constants. These / and d symbols allow us to define both 
antisymmetric and symmetric products among real vectors a, b, • • • in real eight dimensional Euclidean space 1Z S , the 
result in each case being another such vector jlO[ : 

(a A b) r = f rst a s b t , a A b = -b A a; 

(a* h) r = Vid rs tasbt, a*b = b*a. (2-6) 

The significance of these definitions is that aAb and a*b transform just as a and b do, under the eight dimensional 
adjoint representation of SU(3). The matrices of this representation are defined similarly to eq.( pT^ ): 

A G SU(3) -> D rs (A) = ^Tr(X r AX s A^), 

D(A')D{A) = D(A'A), 

D(A) G 50(8). (2.7) 

However in contrast to the SU{2) — 50(3) case, here the matrices D(A) that arise are only an eight-parameter family, 
and so a very small portion indeed of the full twenty-eight-parameter group 50(8). In any case the required properties 
of the products (2.6) are: 

D(A)a A D(A)b = D(A)(a A b), 

D(A)a*D(A)b = D(A)(a*b). (2.8) 
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With this background, we can handle general pure state density matrices for three level systems |10 . Given a 
normalised ip G TL^\ we form the density matrix p = ip^ and expand it in terms of the unit matrix and the X r : 

p = ip^ = \(l+VSn.\). (2.9) 



We then find in place of eq (2.2): 

p 1 " = p 2 = p > 0, Trp = 1 n* = n, n.n = 1, n* n = n (2-10) 

Thus each normalised ray for the three level system corresponds uniquely in a one-to-one manner to a unit vector 
n G 5 7 , the seven dimensional unit sphere in 1Z 8 , which moreover obeys the condition n*n = n. The set of all such 

real unit vectors in 1Z 8 , a subset of S 7 , is the analogue of the Poincare sphere for three-level systems. Since it is in 
fact a very small part of S , we give it a special symbol: 

O = {n G K 8 \n.n = 1, n * n = n} C S 7 C K 8 . (2.11) 

This set O is a connected, simply connected four dimensional region contained in S 7 , and its points correspond 
one-to-one to pure states of a three-level system. It is in fact a representation of the coset space SU(3)/U(2) 
Some interesting geometric properties of O may be mentioned. For two unit vectors ip, tp' G TL^\ we find: 

p = W\ P' = VA^ => Tr(p'p) = W,^ 2 = ^(1 + 2n'.n), 

27T 

< Tr(p'p) < 1 < cos -1 (n'.n) < — . (2.12) 

Thus mutually orthogonal vectors in TiS^ do not lead to antipodal or diametrically opposite points on O, but rather 
to points with a maximum opening angle of ^ radians. Indeed, if n G O, then — n ^ O. If one takes the three 
canonical basis vectors of TL^ as usual, they lead to three distinguished points or "poles" on O: 

(1,0,0) 



(0,1, of 



T 


V3 


1 


— * n.3 = 


= — ' 728 = 


2' 




v/3 


1 


-> n 3 = 


- — ,n 8 = 


2' 


(0,0,1 


T 

-> n 8 = - 


-1, 



(2.13) 



each making an angle of with any other. These properties of O may help one make some sort of mental picture of 
this geometrical object embedded in S 7 . 

The action of SU(3) on vectors in Ti^ leads to adjoint action on O: 

A G 577(3) : V' = M n' = D(A)n. (2.14) 

Thus one has here a (small set of) rigid eight-dimensional orthogonal rotations, which will prove convenient later on. 
Moreover, since O is the coset space SU(3)/U(2), this adjoint action of SU(3) on O is transitive. This will also be 
exploited later. General SO (8) rotations of course do not preserve the region O of S 7 . 

For practical calculations it is convenient to introduce four independent local angle type variables which can be 
used as coordinates over (almost all of) O. Let us write a general unit vetor ip € Ti^ as: 

^ = I 1P2 I , = |^i| 2 + M + m 2 = 1. (2.15) 

W 

Then omitting the part of O corresponding to ^3 = ( this is a two-dimensional region, essentially an S* 2 , see below), 
over the rest of O we introduce 8, <fr, XI1X2 in this way fl3|| : 

(V'l) "02) "03 ) = (overall phase) (e lXl sin 8 cos <f>, e l%2 sin6>sin</>, cos0) 

O<0<§, O<0<§, 0<xi,X2<2tt. (2.16) 
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The limits on 9, (f> reflect the nonvanishing of -03, and the fact that the real three-dimensional unit vector IV^I, 1^3 1) 
has non-negative components. Thus 9, <fi denotes a point on the first octant of an S . Given > 0, xi is the phase 
of -01 relative to ^3 ( an d is well defined except when <p = ^); and X2 is the phase of ■02 relative to ips ( and is well 
defined except when (f> — 0). All this is shown in Fig.l. We need to remember that 9,<p determine the magnitudes 
of the components of ip, while xi>X2 give their relative phases. All four taken together determine one point in the 
portion of O with ipa ^ 0. 




FIG. 1. A pictorial description of almost all of the space O (the space of all states of the three-level system). The local 
coordinates 6, <f>, Xi, Xi are such that 9 and (f> define a point in the positive octant of S 2 i.e. < 6, <f> < ^ and for a given point 
on this octant we have a torus defined by two angle variable < xi> X2 < 2-7T. 



We can easily obtain the expressions for n r in these local coordinates. Combining eqs. fl2.9| , 2.16) we get: 
n r = ^-ip^Xrip : 

n = V3 ^ sin 2 sin </> cos if) cos (%2 — Xi) , sin 2 9 sin <f> cos <f> sin (x2 — Xi ) , 

— sin 2 9(cos 2 <p — sin 2 <p) , sin 9 cos 9 cos <f> cos xi , 

— sin 9 cos 9 cos <f> sin xi > sin 9 cos 9 sin <f> cos X2 > 



sin o cos u sm <p sin X2 



1 



(1 - 3 cos 2 



(2.17) 



The recovery of the Poincare sphere 5* 2 for a two-dimensional subspace of is straightforward. Consider as an 
example vectors ip & T~t^ with vanishing third component (just the points of O omitted in the local coordinatisa- 
tion (pl6|) of O): 




cost 
sin ( 



O<0<-,O<xi, X2<2tt 
Then the eight-vector n has only four nonvanishing components: 



(2.18) 
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(m, n 2 , n 3 ) = — (sin 20 cos (x2 - Xi), sin 20sin (x2 - Xi), c os2</>) , 
1 

«4 = ns = ne = nj = 0. (2-19) 

As </> and {\2 — Xi) vary in the appropriate ranges, we see that we obtain a certain sphere S 2 embedded within O, 
centred on the point (0, 0, 0, 0, 0, 0, 0, of radius Xj^-, and contained entirely within the 1 — 2 — 3 — 8 subspace of 



TZ 8 . If we consider two dimensional subspaces in Tt^ different from (2.18), we clearly obtain SU(3) transforms of the 
above situation. All these various S s are off-centre in TZ 8 : indeed their centres lie on a sphere in TZ 8 centred at the 
origin of TZ 8 and of radius ^ 



III. BACKGROUND TO THE GEOMETRIC PHASE AND PANCHARATNAM'S FORMULA 

Consider a general quantum mechanical system whose pure states are described by unit vectors in a complex Hilbert 
space TL of any dimension. The corresponding ray space will be denoted by TZ. Let C be a continuous piecewise smooth 
parametrised curve of unit vectors in TL: 

C = Ms)\ Sl <s<s 2 }cH, (3.1) 

and let C be its image in TZ, likewise continuous and piecewise smooth: 

C = {p(s) = ipis^is^lsx < s < s 2 } C TZ. (3.2) 

In case ip(s2) and ?p(si) determine the same ray, and in particular tp( s 2) — ^(si) in which case C is closed, the image 
C is closed; however in general we need not assume this. The geometric phase associated with C is the difference 
between a total (or Pancharatnam) phase and a dynamical phase, each of which is a functional of C jBj: 

ip g [C\ = <p p [C] - ipdyn[C], 
ip p [C] = arg(V>(si), V(sa)), 

</?dyn[C]=Im J ds(tp(s),ip(s)). (3.3) 

Si 

The quantity <p g [C] is invariant under both local smooth phase changes in ip(s), and under smooth reparametrisations 
- for these reasons it is a geometric quantity dependent on C rather than on C. 

In this context an important role is played by geodesies in the space TZ Given the continuous curve C C TZ( 
with nonorthogonal end points for definiteness), a nondegenerate positive definite length functional C[C] can be set 
up, which is also reparametrisation invariant: 

S2 i 
C[C) = J ds[(i:(s)J(s))-^(s)J(s))^( S )^(s))Y . (3.4) 

Sl 

Extremising this functional ( with fixed end points), we arrive at the concept of geodesies in TZ. Any Hilbert space 
lift of such a geodesic, with any choice of parametrisation, may then be called a geodesic in TL. It then turns out that 
every geodesic in TZ has vanishing geometric phase: 

C — geodesic in TZ =>■ f g [C] = 0, (3.5) 

and this accounts for their importance. The simplest description of a geodesic, which can always be achieved, is as 
follows. Let the end points of C be and p^ 2 \ assumed nonorthogonal, and choose unit vectors ip^^tpW such 
that: 

(-0 (1) ,-0 (2) ) = real positive. (3.6) 



G 



Thus ipt 1 ' and are "in phase" in the Pancharatnam sense. Then the geodesic C geo C 72. connecting p^ to is 
the ray space image of the following curve C gco C Ti'- 

Cgco = {V>(s)|0 < s < s } , 
V'(s) = "0(0) cos s + 0(0) sin s, 

0(0) = 0«, 0(0) = (> 2) - ^ (1) (0 (1) > (2) )) / (l - (0 (1) ,0 (2) ) 2 ) ' , 

s o = cos- 1 (0«,V (2) )- (3-7) 

Exploiting the fundamental result ( |3.5| ) we obtain a very attractive expression for the geometric phase in the 
following particular situation. Choose a set of points p^ x \ p^ 2 \ ■ ■ ■ , p' n J £ 72 in a definite sequence, assume for 
dcfinitcness that no two consecutive points are mutually orthogonal, and also that p^ and p^ are nonorthogonal. 
Connect p^ 1 to p^ 2 \ p^ to p^ 3 \ • • • , p^ n ' to p^ by geodesic arcs, so that we obtain a closed curve C c 72. in the form 
of an n-sided polygon made up of geodesic pieces. Then we have 0: 

C = geodesic polygon in 72. with vertices p^\ p^ 2 \ • • • , p^") : 
^[C] = -arg(^ (1) ,V (2) )(V' (2) ,V' (3) )---(V' ( " ) ^ (1) ) 
= -argTr(p (1 V (2) ---p (n) ), 
p« = ^W^(i)t, p( 2 ) = V 2)t , • • • , = ^<»ty(»)t. (3.8) 



Here it is evident that the phases of the vectors tp^ ,tp^ 2 \- ■ ■ ,ip^ can be freely chosen. The result ( |3.§| ) con- 
nects the geometric phase for a closed polygon to the Bargmann invariant of quantum mechanics, the expression 
(^( 1 ),V> (2) )(0 (2 \0 (3) )---(0 ( ™ ) ,0 (1) ), namely: the for mer is the negative of the argument of the latter. The point to 
emphasize is that the definition of the Bargmann invariant requires specifying just the vertices of the polygon, while 
the definition of the geometric phase requires also connecting them in sequence by geodesic arcs so that we have a 
closed loop C c 72. 

With this background from the general theory of the geometric phase we relate these results to the case of two 
level systems, quote the Pancharatnam formula and then give the general expression for <p g [C] for three-level systems. 

For two-level systems we have seen that the space 72 is the Poincare sphere S . It is now a happy coincidence that 
geodesies in 72 map exactly on to geodesies on S 2 in the more familiar Euclidean sense. This can be seen as follows. 
Without loss of generality, by using a suitable SU(2) transformation, we may assume that the points p^\ p' 2 ' to be 
connected by a geodesic are the points n^ 1 ) = (0, 0, l),n^ 2 ^ = (sin2a, 0, cos 2a) on S 2 , with representative vectors 
^(1)^(2) e 7^(2) chosen as follows: 

Then applying the result (3.7), the ray space geodesic connecting p^> to p^ is determined as follows: 



V>(s) 



cos s 
sin s 



n(s) =Trp(s)a= (0>(sW(s)) 

= (sin 2s, 0, cos 2s) , < s < a. (3.10) 

We see that the curve described by n(s) on S 2 is indeed a great circle arc, a part of the "Greenwich meridian"; and 
by the action of SU(2) on 72 translated into the action of 5*0(3) on S 2 , we conclude that any ray space geodesic in 
72 appears as some great circle arc on S 2 . Thus the two definitions of geodesies do coincide in this case. 

Let now A, B, and C be any three points on S 2 , no two being diametrically opposite to one another. Joining them 
by great circle arcs (each less than 7r in extent) we get a geodesic triangle A (A, B, C) on S 2 . Then the Pancharatnam 
formula || Q Q is the statement that for any two level system: 

^[A(A,B,C)] = iQ ) 

il = solid angle subtended by the triangle A, B, C at the origin of S 2 . (3-11) 

Here the right hand side is interpreted to be positive(negative) if, as viewed from outside of S 2 , the triangle ABC is 
described in the counter clockwise (clockwise) sense. It is this formula that we shall generalise in Section V. 
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Now we give the general formula for geometric phases for three level systems |Hj. Consider a closed loop C <Z O, 
and assume for defmiteness that ^3 7^ throughout. Then <^ g [C] is given by the following integral along C: 



tp g [C] = - f sin^6»(cos 2 (/)dxi +sin 2 (/)dx2) (3.12) 
ceo 

We note that this formula holds when C is a closed loop. In the succeeding Sections we develop the properties of 



geodesies in O and then generalise eq. (3.11) 



IV. RAY SPACE GEODESICS FOR THREE-LEVEL SYSTEMS 



We have seen that the transitive action of ST/ (3) on the ray space O for three- level systems is given by (an eight- 
parameter subset of) rigid orthog onal rotations in real Euclidean eight dimensional space, when points of O are 



identified with vectors n as in eq (2.11). Based on this we may describe the details of any one conveniently chosen 
geodesic in O; and then any other would be a suitable SO(8) transform of this one, so that the geometrical shape 
and structure in an intrinsic sense are unaltered. 

Guided by the constructions of Section III, eq. (3J3, 3.10| ) let us choose two points in O correponding to the following 
two unit vectors in TL^: 

n (1) = (0,0,0,0,0,0,0 
n (2) =^^0,0, - sin 2 a, 0,0, 2 sin a cos a, 0, i(l - 3cos 2 a)^ (4.1) 

It is easy to see that, given any pair of three level pure state density matrices and p^ such that Tr(p^ p^) = 
cos 2 a > 0, we can exploit the action of SU(3) on H^ 3 ', and freedom of phases, to put and p^ into the 
configurations corresponding to the vectors ip^ 1 " 1 and ip^ a bove . Then the ray space geodesic C g eo connecting 
and n^ 2 ) is easily found on the basis of the general formula (3.7): 





C7(°>: n( S ) = ^( s )t A ^( s ) 

= ^-(0,0, -sin 2 s ,0,0, 2sinscoss, 0, -=(l - 3cos 2 s) ) , < s < a. (4.2) 
2 V3 / 

This curve {n(s)} C O has only three nonvanishing components, namely, n3(s),ne(s) and ris(s). The interesting 
questions are whether it is a plane curve, whether it coincides with a geodesic arc as defined in the sense of eight 
dimensional Euclidean geometry on S 7 , and whether it has any other important geometric features. 

The first observation we may make is that since the components ni, 712, 714, n^, all vanish, this curve Cgeo lies 
entirely in a three dimensional subspace of 1Z S . By SU (3) action this statement is then true for all ray space geodesies 
C gco when drawn in O. Next we remark that a geodesic on S 7 in the Euclidean sense would be part of the intersection 

of a two dimensional plane in 1Z 8 passing through the origin, with S 7 . We can immediately see that C g c is not of this 
kind; for example, the three vectors n(0), n( ^), n(a) are easily checked to be linearly independent. Thus in striking 

contrast to the situation on the Poincare sphere S 2 for two level systems, here ray space geodesies on O are not 
geodesic arcs in the sense of Euclidean eight dimensional geometry. 

Nevertheless we can show that C g o is a plane curve, though the plane on which it lies is off center in 1Z S , i.e. it does 

not pass through the origin or the center of S 7 . For this we need to work with the combinations ^Zlzi^+JM ; _3 — i \ // 3n$ 
obtained from n^,ns by an orthogonal transformation. Then we find: 
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V3n 3 (s) + n s (s) _ 1 
2 ~ ~2 ' 

'^(s) = VSsin s cos s , 

"3W~ 2 >/3n 8 W = _^ (1 _ 2cos 2 8) 



(4.3) 



Thus, in the three dimensional 713 — ri6 — subspace, this geodesic C g cl is a curve lying in the two dimensional plane 

= -|. Ray space geodesies C gco for three level systems, when described as curves { n(s)} C O, are really 
like arcs of constant latitude on S 2 , and not geodesic arcs at all. For this reason, hereafter the term geodesic will refer 
exclusively to ray space geodesies, with no further qualification. 

In the description of O using coordinates 9, <f>, Xi> X2 the geodesic C S eo appears as follows: 

:6{s) = s, = |, xi(s)undefmed, = 0, < s < a. (4.4) 

Thus the projection of C g co onto the octant of S 2 with spherical polar angles 9, <fi happens to be a great circle arc in 
the usual sense. However this is not expected to be true for a general geodesic C gco C O. 

One can now ask for the most gener al th ree level system Hamiltonian which reproduces the evolution of ip(s) with 
respect to s, as already specified in cq. (4.2), and whether it can be independent of s. A general hermitian Hamiltonian 
can be written in terms of nine real functions of s as: 

H(s) = h (s) + h(s).A, (4.5) 

while the quantum mechanical evolution equations for tp(s) and p(s) become: 

±n(a) = 2h( S ) A n( S ). (4.6) 
as 

An easy calculation shows that the most general H(s) has four independent real arbitrary functions of s, namely 
a(s), b(s), c(s), d(s) in: 

h(s) = (a(s) cos s, b(s) cos s, v / 3c(s) + <i(s)(cos 2 s — sin 2 s), 
— a(s) sins, — b(s) sins, d(s) cosssins, —1, c(s)) 

2 

h (s) = —j=c(s) - d(s) sin 2 s. (4.7) 
v3 

Moreover we also find with this general Hamiltonian : 

Tr(p(8)H(8)) = M8),H(8)1>(8)) 

2 

= h o(s) + -^=n(s).h(s) 

= 0. (4.8) 



This is consistent with the vanishing of the dynamical phase as the state evolves along {ip(s)} of eq ( J4.2| ), which is 
directly checked to be true. If we in addition wish to have a constant Hamiltonian producing this evolution, we must 
set a(s) = b(s) = d(s) = 0, c(s) = cq = constant and then we get: 

#c ( o°L. = (4=1 + V3A 3 + A 8 ) c - A 7 . (4.9) 

It is possible to express this simple constant Hamiltonian directly in terms of the end points n^, n' 2 ) of the geodesic 
Cgco given in eq. (4.1). Namely one finds: 

(n (1) A n (2) ) r = -~ (sinacosa)^ (4.10) 







which, if we set the constant Co to zero, leads to: 

(0) nf) A nW.A 

"-const. | n (X) An (2)| ' 

We can now generalise these results to any geodesic C goo C O connecting any two points n,n' E O. Writing 
n.n' = i(3cos 2 a — 1), we have the result that the constant Hamiltonian 

Hconst. = 7r~- ' (4-12) 

3 sin a cos a 

can produce evolution along the geodesic C gco ■ This generalises well known results in the case of two- level systems [Jrl . 



V. THE PANCHARATNAM GEOMETRIC PHASE FORMULA FOR THREE-LEVEL SYSTEMS 



The purpose of this Section is to obtain the generalisation of the result ( |3.11 ) for three-level systems. For this 



we must determine the geometrical description of the most general geodesic triangle on O, upto an overall SU(3) 



transformation, and then use the general connection (3.8) to find the geometric phase associated with this triangle. 

Let A, B and C be three general points in O such that no two of them enclose an angle of ^ radians. Using the 
freedo m o f common SU(3) action we can transport A to a position A^ with a representative vector of the form ip^ 



of eq (|4.l| ). In this process let B and C move to locations B' and C: 

SU(3) , . 
A,B,C^ A(°\B',C 

A<® - p<U = V (1) V (1)t , = | | (5.1) 

This p^'(ip^) is invariant under an U(2)(SU(2)) subgroup of SU(3). Exploiting this we can next transport B' to a 
position while leaving A^ fixed, and simultaneously C to some C" , such that: 

AW,B',C ,V S J ] AW,B«>\C" : 
£?(°) p W = ^V 2)t , ^ 2 )= | sin^ ,0<e<|, (5.2) 

thereby introducing the angle £. We have also secured that (tp^,^ 2 ') is real positive. Now these two density matrices 
pv 1 ), p( 2 ) are invariant under a particular diagonal U(l) subgroup of SU(3), whose elements are: 

d(/3) = diag (e~ 2l/3 , e i/3 , e i/3 ) € SU(3), < (3 < 2tt. (5.3) 

We now have the freedom of transformations e?(/3), which leave A^ and B^ unchanged, to move C" to a convenient 
position . A little thought shows that this can be achieved as follows: 

AW,B®\C'™A<®,B®\C<® : 
sin i] cos £ \ 

C(o) ^(3) = [ e i X2 sinrjsinC , < r\ < ~, < C < |, < X 2 < 2tt. (5.4) 
cos?? / ~ 




We have parametrised tp^ in the manner of eq. ( 2.1 6| ) : the U(l) freedom allows us to transform the phase \i to 



zero, and an overall phase freedom has been used to make (^W,^ 3 )) real positive. We now see that the description 
of the most general geodesic triangle in 0, upto an overall SU(3) transformation, involves the four angle parameters 
?i 7 7jCiX2- These are therefore intrinsic to the shape and size of the triangle. This counting agrees with the fact 
that SU(3) is eight dimensional, and choosing three points on O independently involves choosing twelve independent 
coordinates. 
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Now we consider the geometric phase for the geodesic triangle A' ', B^°\ . Us i ng t he result (3.8) based on the 
Bargmann invariant, and the conveniently chosen representative vectors in eqs. ( |5.l| , 5.2 , 5.4), we find: 

= +arg(V' (3) ,V' (2) ) 

= arg(cos £ cos r/ + sin £ sin r\ sin £e _IX2 ) 

= arg(l + tan^tan?7sinCe _lX2 ) (5.5) 



This is the generalisation of the Pancharatnam formula (3.11). We see that the phase xi plays an important role: 
(p g [A^B^C^] can be nonzero only if \2 is nonzero. 

We can also express this geometric phase directly in terms of the vectors n^ 2 \ S O corresponding to the 
vertices A(°), C(°) of the geodesic triangle: 



<p g [A^B^C^} = -arg Tr(p^p^p^) 



= — tan 



2V3n( 1 ).n( 2 ) A n( 3 ) 



(nW + n( 2 ) + n( 3 )) 2 + 2n( 1 ).n( 2 ) * n( 3 ) - 2 



(5.6) 



In this form the SU(3) invariance is explicit. 

We may collect our results of this and the previous Section to say: given any geodesic triangle in the ray space O 
for three-level systems, it is possible to find a piecewise constant Hamiltonian to produce evolution along this triangle, 
such that the dynamical phase vanishes, and the geometric phase is then given by the intrinsic expressions (5.5, (5.6). 



VI. CONCLUDING REMARKS 



We have exploited the newly constructed extension of the Poincare sphere representation from two to three-level 
quantum systems, to develop a generalisation of the Pancharatnam geometric phase formula for cyclic evolution of 
such a system along a geodesic triangle in state space. We have found that such a triangle is intrinsically defined by 

four angle type parameters in contrast to the three needed for defining a triangle on the Poincare sphere S 2 ; and have 
obtained the explicit and simple expression for the geometri c pha se in terms of them 

It is easy to check that the original Pancharatnam formula (3.11) emerg es fo r m e q (5.5) as a particular case. Namely 



if we take C = \i the three Hilbert space vectors tpW ,ip@) ~ipW of eqs (5.1, 5.2 , 5.4]) all lie in the two dimensional 
2 — 3 subspace of Ti^ , and involve just three angle parameters £, r), \2- 



(1) _ 



r r2 ' = \ sin£ J , i/> (3) = ( e** 2 sin 77 
cos£ / y cos Tj 



(6.1) 



Let the corresponding rays be represented by unit vectors ni, 112, 113 on an S 2 and let the sides of the corresponding 
spherical triangle be a, b and c. Then from eq (|2.l|), we have: 



ni.n 2 = cosa = 2|(^ (1) , ^ (2) )| 2 - 1 = cos2£, 
ni.n 3 = cosb = 2|(-0 (1) , V (3) )| 2 - 1 = cos2?7, 

n 2 .n 3 = cose = 2|(V> (2) , "0 (3) )| 2 - 1 = 2| cos£cos?/ + sin £ sin ?7e lX2 1 2 - 1. 



Thus a = 2£, b = 2r\ and in the case of the angle c we have 

c 



cos — = I cos £ cos 77 + sin £ sin r/e 1 * 2 \ 



Now from eq (5.5) the geometric phase in this case is given by: 

<p g [A^B^C^] = ar ff (cos£cos77 + shi<esin77e- iX2 ), 
or equally well (apart from the sign) by: 



(6.2) 



(6.3) 



(6.4) 
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o«t<p g [A m B<f»cM]) = gggj£gg^+j£i^g£gg2Q, (6 . 5) 
y |cos^cos77 + sin^sin?7e l X 2 | 

Focussing on the \2 dependence here (and as we have seen earlier this is the crucial aspect), and using eq. ( |3.2| ) we 
have: 

... (1 + cos a + cos b + cose) . „. 

cos £ cos 77 + sin £ sin r? cos % 2 = -. 77. tjt, . (6-6 

4 cos a/2 cos 6/2 



leading to 



iM A {0) B^C^}) = (l+^ + ^ + cosc) ^ 

V ^ 9L v 4cosa/2cos&/2cosc/2 v ; 



However the right hand side here is precisely cos |A(a, b, c), where A(a, 6, c) is the solid angle subtended at the origin of 
S 2 by the spherical triangle with sides a, b and c. In this way the Pancharatnam result (p g [A^B^C^] — |A(a, 6, c) 
is recovered. (The sign can also be recovered with some additional effort). 

This same verification leads us to the following significant remark: as long as one is interested in cyclic evolution of 
any quantum s yste m along a geodesic triangle in ray space, however large the dimension of the Hilbert space TL may 



be, our result (5.5) for the geometric phase is applicable and is completely general. This is because a triangle involves 
(at most) three independent vectors vj^\ ip^ , £ Tt, and these always lie in some three dimensional subspace 
of 7i. Thus nothing additional is needed to handle geometric phases for cyclic evolutions along geodesic triangles 



for iV-level systems, for any N > 4. In this sense the original Pancharatnam result (3.11) deals with a degenerate 



case, since it is concerned with evolution along a triangle but for a two-level system. There is then no place for the 



additional angle variable £ to appear as it does in the fully general formula (5.5). 

A discussion of feasible experimental schemes to check the validity of our main result (5.5), at least in some nontrivial 
cases which do go beyond the two-level situation, will be the subject of a forthcoming publication. 
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